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Abstract 

We propose a new spectral model for the Seyfert 1 Galaxy MCG-6-30-15 that can describe most 
of the 1-40 keV spectral variation effectively by a single parameter. Our spectral model includes 
three continuum components; (1) a direct power-law component, (2) a heavily absorbed power-law 
Cln component by thick photoionizcd material, and (3) a cold disk reflection component far from the 

black hole with moderate solid-angle (J7/27T « 0.3) including a narrow fluorescent iron line. The 
first two components are affected by two warm absorbers with different ionization states. The 
ionized iron K-cdge of the heavily absorbed power-law component can explain most of the seemingly 
broad "disk line" spectral feature, so that our model does not require an additional broad iron 
emission line. The observed spectral variation is fully described by change of the three spectral 
parameters; normalization of the direct power-law component, that of the heavily absorbed power- 
law component, and column density of the lower-ionized warm absorber. Furthermore, normalization 
of the direct power-law component is strongly anti-correlated with that of the heavily absorbed 
power-law component, as well as with column-density of the lower-ionized warm absorber; the two 
constraints among three free parameters suggests presence of a single principle parameter. In fact, in 
the course of observations, the "total normalization" , which is sum of the normalizations of the direct 
power-law component and the absorbed component, is not significantly variable, while the "covering 
fraction" , which is defined as ratio of the normalization of the absorbed component to the total 
normalization, is much more significantly variable. We propose that MCG-6-30-15 has a moderately 
extended, not significantly variable central X-ray source, partially covered by intervening optically 
thick clouds having internal ionization structure, so that observed flux and spectral variations are 
primarily caused by variation of the partial covering fraction. These absorbing clouds probably 
correspond to the fast-moving broad line region (BLR) clouds. 
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1. Introduction 
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The Seyfert 1 galaxy MCG-6-30-15 is the primary "disk-line" target, in which a relativistically broadened 
iron K-emission line is suggested to originate in the innermost region of the accretion disk around the black hole 
(e.g., Tanaka et al. 1995). MCG-6-30-15 is also known to exhibit characteristic X-ray spectral variations; Inoue 
& Matsumoto (2001), Fabian et al. (2002) and Matsumoto et al. (2003) reported significantly small variability 
in the iron line energy band of MCG 6-30-15. It is found that the energy-dependent Root Mean Square (RMS) 
variability is suppressed around the iron line energy, particularly in longer timescales. Inoue and Matsumoto 
(2001, 2003) proposed that the absorbed spectrum due to photo-ionized warm absorbers mimic shape of the 
strongly red-shifted iron line, and variation of the warm absorbers may explain the apparent small variability 
in the iron line energy band. 

On the other hand, Miniutti and Fabian (2004) proposed that the suppressed variability in the iron line 
band may be explained by the general relativistic light-bending effects, which take place in the very vicinity 
of the black hole. Miniutti et al. (2007) claim that the Suzaku spectrum of MCG-6-30-15 is in fact consistent 
with the light-bending model. Niedzwiecki & Zycki (2008) and Niedzwiecki & Miyakawa (2010) independently 
re-examined the light-bending model, and concluded that it is possible to explain the suppressed variability in 
the iron energy band assuming a particular configuration and movement of the illuminating source. Zycki et al. 
(2010), on the other hand, claim that the light-bending model is not able to explain the Suzaku spectrum, if a 
broader energy range is adopted than used by Miniutti et al. (2007). 

The primary question in MCG-6-30-15 is whether the iron line is truely broad and its little variability is a 
consequence of the general relativistic light bending effects or not. The critical point to address the problem is 
to model the underlying continuum spectrum correctly, since the broad iron line parameters are dependent on 
the choice of the continuum spectral model. In fact, Miller, Turner, and Reeves (2008) claim that the disk-line 
is not required, when the continuum spectrum is modeled with multiple warm absorbers. 

Recently, using the Suzaku data taken in January 2006, Miyakawa et al. (2009) (hereafter Paper I) found 
a clear correlation between the intensity in the 6-10 keV band and the spectral ratio of 0.5 - 3.0 keV / 
6-10 keV in a time-scale of < 10 5 sec. Essentially the same spectral correlation was confirmed for longer 
timescales up to ^14 years in the RXTE data (Miyakawa 2010). Paper I adopted a spectral model composed 
of a direct power-law component, its reflection component by a neutral matter (with the solid-angle f2/27r« 1), 
two warm absorbers with different ionization states. The observed spectral variation requires change of the 
apparent slope of the direct component, whereas the shape and intensity of the reflection component being 
invariable. A mildly broad iron emission line was required at 6.42 ±0.06 keV with an intrinsic width of 0.29 keV 
(1 u) and the equivalent-width 100 cV; such an extremely broadened "disk line" as is claimed by Miniutti et al. 
(2007) was not required because of difference of the continuum models. In Paper I, we confirmed a well-known 
observational fact of MCG-6-30-15 that the reflection component as well as the iron line flux is hardly variable 
while the direct component is so variable, which has yet to be explained. 

In this paper, following Paper I, we attempt to comprehensively understand spectral variability of MCG-6-30- 
15 using Suzaku and Chandra archival data. Our goal is to find a spectral model which can naturally explain the 
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observed spectral variation of MCG-6-30-15. If we are successful to construct such a reasonable spectral model, 
we will see if the relativistically distorted iron emission line is truely present or not in the energy spectrum. We 
believe this is an important step forward to fully understand X-ray spectral properties of MCG-6-30-15, Seyfert 
1 galaxies, and black holes. 

2. Observation and Data Reduction 

We use the same Suzaku data used in Paper I, taken in 2006 January, between 9-14 (143 ksec exposure), 
23 -26 (99 ksec) and 27 - 30 (97 ksec). For Suzaku data reduction, sec Paper I. 

The Chandra satellite observed MCG-6-30-15 several times. In this paper, we use the High Energy 
Transmission Grating Spectrometer (HETGS) data taken in 2004 May, where the source was observed four 
times between May 19 and 27 resulting in a good exposure time of 522 ksec. For data reduction we used 
the CIAO 4.0 software package; "tgextract" is used to produce PHA2 spectral files from the Level2 event file. 
HETGS is composed of High Energy Grating (HEG) and Medium Energy Grating (MEG). The +1 and -1 
orders of the HEG were combined, and so were the +1 and —1 orders of the MEG. The HEG and MEG spectra 
were fitted separately. 

3. Data Analysis and Results 

3.1. Introduction of the "3- component" model 

In Paper I, spectral variation of MCG-6-30-15 was studied assuming the "two-component" spectral model 
given as, 

F = W H W L (N D + N R R)P + I Fe , (1) 

where P is the intrinsic power-law spectrum with a high-energy cut-off, No is normalization factor of the 
directly observed component, Nr is normalization factor of the reflected component, R is reflection albedo by 
the cold optically thick matter, Wh and Wl represent attenuation by a high- and a low-ionized warm absorbers, 
respectively, and Ip e is an iron K a emission line. Effect of interstellar absorption is always taken into account 
in the spectral fitting but its term is omitted for simplicity in the above expression and hereafter. 

Paper I found that the spectral variation at various timescales is successfully explained by variations of only 
two parameters, normalization factor of the directly observed component (Njj) and apparent slope of the direct 
component, such that the slope gets steeper for larger normalization. The apparent slope change is described by 
either change of the index of the power-law component (P), ionization state of the low- ionized warm absorber, 
or column density of the low-ionized warm absorber (Wl). 

Below, we scrutinize the observed energy spectra more carefully, and introduce a more elaborated model. 
Using the Chandra HETGS, Young et al. (2005) resolved a narrow emission line at 6.4 keV in MCG-6-30-15, 
which we confirmed via re-analysis of the archival data. The equivalent width of the line is measured as ~ 20 
eV, which requires a cold reflector with a solid angle of Q/2tt ~ 0.3 viewed from the central X-ray source (e.g., 
George & Fabian 1991). On the other hand, the best-fit values of Nl> and Nr obtained in Paper I indicates that 
Nr/Nd = £l/2ir sa 1. Consequently, in order to account for the observed disk reflection spectrum of an amount 
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of S1/27T ~ 1 including the outer disk reflection with fi/27r ~ 0.3, we need another disk reflection component 
corresponding to an amount of 0/27T ~ 0.7. 

Thus, we separate the Nr term in (1) into two terms. One is the reflection continuum component from a 
cold reflector emitting a narrow iron K a line, solid angle of which is D,/2ir ~ 0.3. The other is introduced to 
represent the remaining part of the Nr term but is generalized to cope with several possibilities: First, this 
reflector might exist in the vicinity of the central X-ray source and could be ionized. Second, the reflected 
spectrum looks similar to an absorbed spectrum and thus the seemingly reflection term might be due to in 
reality partial absorption (e.g., Matsuoka et al. 1991). Third, this term could be variable independently of other 
components. 

Resultantly, we introduce the new model expressed as 

F = W H W L (N 1 +N 2 W 2 )P + N 3 RP + I Fe , (2) 

where Ni is normalization of the direct component, N 2 is for the reflection component or absorbed component, 
W 2 stands for another warm absorber which approximately represents energy dependent absorption or albedo 
of the ionized optically thick matter. Normalization factor of the outer reflection component, N 3 is constant 
and related to the average of Ni so that N3/ < Ni >= H/2tt ~ 0.3, and Ip e is a narrow iron emission line with 
a fixed equivalent width of 20 eV. Hereafter, this model is referred as the "3-component" model. 

We try the 3-componcnt model to the Suzaku average spectrum. We found the fit improves if we additionally 
put a weak iron neutral edge at 7.11 keV and an absorption line at 7.0 keV. The former might be due to a 
slight difference of the assumed f2 from the exact value, and the latter due to deficiency of the H-like iron 
absorption line equivalent-width already included in the high-ionized warm absorber (calculated with XSTAR; 
see Paper I). Figure 1 shows the fitting result for the Suzaku average spectrum with the 3-component model 
in 1 - 40 keV. Parameters of the high-ionized and low-ionized warm absorbers are similar to those in Paper I 
(Table 1). For the newly introduced reflection or absorption component (W 2 ), we obtain Nr ~ 1.6 x 10 24 cm~ 2 
and log£ = 1.57. For the time being, we leave its interpretation as either disk reflection component or heavily 
absorbed component. Ratio of the normalization of the heavily absorbed component to the direct component is 
0.60/1.44 w 0.4. If interpreted as disk reflection, this suggests Vl/2ii ~ 0.4. Or, in the absorption interpretation, 
the source is partially covered by the ionized thick matter by 0.6/(1.44 + 0.6) «30 %. 

3.2. Examination of presence of the "disk line" 

Several authors claim presence of an extremely broad "disk line" in the energy spectrum of MCG-6-30-15 
(e.g., Miniutti et al. 2007). Therefore, we examine if presence of the disk line is reconciled with our 3-componcnt 
model. We analyze the time-average spectrum of MCG-6-30-15 adding a disk line model from a fast-rotating 
Kerr black hole (Laor 1991). Thus, the model we adopt here has both the narrow line from very far from the 
black hole, and a putative broad line which is from the region very close to the black hole. We try to fit the 
spectrum with Laor's disk line model with varying the inner disk radius, r !r j. 

Table 2 shows the best-fitting result for the Suzaku average spectrum adding Laor's model to the 3-component 
model, where the fit is slightly improved with ri n 200 r g and the reduced chi-squarc is 1.07 (x 2 /d.o.f = 
214/200). Note that the additional line is only mildly broadened with r,„ w 200 r g , and the equivalent width 
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Table 1. Results of spectral fitting in 1-40 keV for the Suzaku XIS/PIN average spectrum with the 3-component model. 



Interstellar absorption 




N H (10 21 cm- 2 ) 


1 -°-0.2 


W H 




N H (10 23 cm- 2 ) 


2 4+ 1 - 8 


logC 


3.37±0.04 


W L 




N H (10 21 cm" 2 ) 


3.7±1.3 


log? 


1 K4+0.25 


JVi 




10" 2 ph/s/cm 2 at 1 keV 


1.44±0.03 


N 2 




10- 2 ph/s/cm 2 at 1 keV 


0.60±0.09 


W 2 




N H (10 24 cm- 2 ) 


1 fi +0 - 6 

J -°-0.1 


logC 


i r 7 +0.21 

-0.19 


P 




photon index 


1 01 +0.02 


E cut (keV) 


160(fixed) 


N 3 


7Vix0.3 


i?t 




coslncl 


0.866(fixed) 


iFe 




line E (keV) 


6.35 (fixed) 


sigma (keV) 


0.01 (fixed) 


norm (10 _5 ph/s/cm 2 ) 


1.2±0.3 


EW (cV) 


27±7 


edge E (keV) 


7.11 (fixed) 


MaxTau 


0.05±0.02 


line E (keV) 


7.0 (fixed) 


sigma (keV) 


O.Ol(fixed) 


norm (10 _6 ph/s/cm 2 ) 


-6.8±2.4 


line E (keV) 


2.35±0.02 


sigma (keV) 


O.Ol(fixed) 


norm (10 _5 ph/s/cm 2 ) 


-2.3±0.5 


reduced chi-square (d.o.f) 


1.12 (203) 



f "pexrav" model in xspec is used. 
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Fig. 1. Spectral fit result for the 1-40 keV average spectrum for Suzaku XIS/PIN with the 3-component model. 
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Fig. 2. Chi-squares of fitting to the average spectrum including Laor's disk line model, shown as a function of the inner 
radius. 
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Table 2. Results of spectral fitting in 1-40 keV for the average spectrum adding Laor's disk line to the 3-componcnt model. 



Interstellar absorption 




N H (10 21 cm- 2 ) 


O r+1.4 

°-°-0.5 


Laor 's disk-line 




E (kcV) 


6.40±0.06 


Index 


-3 (fixed) 


n„ (Rg) 


22n +iso 


Tout (Rg) 


400 (fixed) 


inclination (deg) 


30 (fixed) 


norm (10~ 5 ) 


1 O+0.4 

^-o.e 


EW (cV) 


40±20 


W H 




N h (10 23 cm" 2 ) 


o n+0.9 
z - J -0.4 


logC 


3.35±0.04 


W L 




N H (10 21 cm" 2 ) 


3 7+1-4 

°- ' -0.8 


logC 


-, r 9 +0.21 
1 - oz -0.12 


Ni 




10- 2 ph/s/cm 2 at 1 keV 


i '^^-0.02 


N 2 




10" 2 ph/s/cm 2 at 1 keV 


^+ 05 


W 2 




N H (10 24 cm" 2 ) 


1 -' -0.2 


log £ 


1 - DO -0.16 


P 




photon index 


1 OT+0-02 
1 - yi -0.01 


E cut (keV) 


160(fixed) 




iVixO.3 


R) 




coslncl 


0.866 (fixed) 






line E (kcV) 


6.35 (fixed) 


sigma (keV) 


0.01 (fixed) 


norm (10 -6 ph/s/cm 2 ) 


K 4+2-4 


EW (eV) 


iz -10 


edge E (keV) 


7.11 (fixed) 


MaxTau 


0.04±0.02 


line E (keV) 


7.0 (fixed) 


sigma (keV) 


O.Ol(fixed) 


norm (l()- 6 ph/s/cm 2 ) 


-6.1+2S 


line E (keV) 


2.35±0.02 


sigma (keV) 


O.Ol(fixed) 
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of the broad disk line is just 40 ± 20 eV. Figure 2 shows the chi-squares of fitting to the average spectrum as a 
function of the inner disk radius. As seen from this figure, A chi-squares become greater as rj n gets less than ss 
200 r gi and it exceeds 2.71 when ri n is less than 9 r g . Namely, our 3-component model is only reconciled with 
a mildly broad disk line from the region r in >9r 9 with 90 % confidence. 

The reason we do not require a strong, extremely distorted disk line is primarily that the newly introduced 
heavily absorbed component in our 3-component model (the N2W2P term in Equation 2) has an ionized iron 
edge feature which resembles the disk line shape (Figure 1), and secondly that the warm absorbers (Wl and Wh 
in Equation 2), which are already included in Paper I, have spectral curvatures below ~5 keV without which 
the residual may look like a low-energy tail of the disk line. Consequently, we see that the disk line feature is 
dependent on the choice of the continuum spectral models, and an extremely distorted disk line is not required 
in our 3-componcnt model. 

3.3. Spectral fitting for the sliced spectra 

The spectral model fitting alone does not tell if the heavily absorbed component we introduced in the 3- 
component model, the N2W2P term in Equation 2, is really an absorbed component or due to disk reflection. 
To answer this question, we further study spectral variations observed with Suzaku, using the 3-componcnt 
model as a working hypothesis. 

We apply the 3-component model to the eight intensity-sliced spectra made in Paper I. Here, all the param- 
eters are the same for the eight spectra except the following three parameters; the direct power-law normaliza- 
tion, N±, the absorbed power-law normalization, N2, and ionization degree of the low- ionized warm absorber, 
Wl- Attenuation by a warm absorber is expressed as Wl ~ exp(— a(E,^) Nh,l), where cr(E,£) is energy and 
ionization-degree dependent cross-section and Nh.l is the hydrogen column density of the low-ionized absorber. 
Its variation is due to cither a variation of a(E,£) or Nh,l, or both. Here, we investigate two cases, one with a 
variation of the cross-section a(E,^) due to change of the ionization parameter £, the other with a variation of 
N h ,l- 

Eight spectra are fitted simultaneously, and the both cases equally exhibit acceptable fits with reduced chi- 
squarc values 1.09 and 0.95 (d.o.f = 1135), respectively. Thus, we could not tell if the ionization parameter or 
column density of the low-ionized warm absorber is more variable from spectral fitting alone. 

When Nh.l is fixed, the top and bottom panel of Figure 3 shows the relation between N\ and log£ and that 
between N\ and N2, respectively. Similarly, when £ is fixed, the top and bottom panel of Figure 4 shows the 
relation between N± and Nh,l and that between N\ and N2, respectively. In either case, correlations between 
each two of the three parameters are obvious. Note that we sec two independent correlations among three free 
parameters. This indicates that there should be a single principle parameter which is primarily responsible for 
the observed spectral variations (section 4.1). 

3-4- Spectral fitting for Chandra/ RET GS data 

We have seen that the 3-componcnt model is effective to describe the Suzaku spectral variation, and that 
the ionization degree of the low-ionized warm absorber is variable with intensities. Next, we will see if the 
3-component model is valid to describe the Chandra HETG spectra and their variations. Strength of Chandra 
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Fig. 3. Correlation among the three parameters, the direct power-law normalization (N±), ionization parameter of the 
low-ionized warm absorber (£), and the absorbed power-law normalization (N2), when applying the 3-component model to 
the eight intensity-sliced spectra. (Top) Relation between N± and £. (Bottom) Relation between Ni and N2. 
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Fig. 4. Correlation among the three parameters, the direct power-law normalization (N\), column-density of the low-ionized 
warm absorber {Njj £), and the absorbed power-law normalization (N2), when applying the 3-component model to the eight 
intensity-sliced spectra. (Top) Relation between JVi and Nh,l- (Bottom) Relation between iVi and N2. 
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HETG compared to Suzaku is that it can study variation of individual low-cncrgy absorption lines corresponding 
to change of the ionization degree. 

First, we fit the Chandra/HETGS spectra with the 3-component model in 1.0 - 7.5 keV. Since the energy 
range is much smaller than that of Suzaku, some parameters are not constrained and thus fixed to the Suzaku 
best-fit values. In Figure 5 and Table 3, we show the fitting result for the average spectrum. We see that the 
Chandra HETG spectra are successfully modeled with the 3-component with similar parameters of those of 
Suzaku. 

Next, we study possible ionization degree variations in the Chandra/HETGS data by extracting the "bright 
spectrum" and "faint spectrum" as follows: 1) Create a light curve with a time-bin-width of 128 sec, and 
calculate average counting rates. 2) Create the "bright spectrum" from the period when the MEG and HEG 
count rates (0.4-10 keV) are higher than the average, and the "faint spectrum" when the count rates are lower 
than the average. 

We fit the bright/faint spectra in the 1.15-1.55 keV band with a single power-law and two negative gaussians 
for the MgXI (He-like 1.34 keV) and MgXII (H-like 1.47 keV) absorption lines. We also fit the bright/faint 
spectra in the 1.75-2.10 keV band separately with a single power-law and two negative gaussians to investigate 
for the SiXIII (He-like 1.85-1.86 keV) and SiXIV (H-like 2.00 keV) lines. In Figure 6, we show the fitting result 
for the bright/faint spectra in both energy bands. We could fit the bright and faint spectra with the common 
photon index, while only the power-law normalization and the normalization of absorption lines are varied. The 
best fit parameters are shown in Table 4. 

We notice that the MgXI absorption line equivalent width is larger in the faint state, while that of the MgXII 
line is larger in the bright state. This is understood as due to change of the ionization state of the warm 
absorber, such that MgXI is more abundant in the faint state, while MgXII is more abundant in the bright state. 
In fact, ion fraction of MgXI and that of MgXII become equal at log£ ~ 2 (Kallman & Bautista 2001), that is 
about the best-fit value of the low-ionized warm absorber (Tables 3.4). 

As for Si absorption lines, the SiXIII line equivalent width is oppositely larger for the bright spectrum, which 
is explained by the low- ionized warm absorber with log£ ~ 2, since ion fraction of SiXIII increases when log£ 
varies from ~ 1.5 to ^2 (see Kallman & Bautista 2001). Fraction of SiXIV should increase more dramatically 
when log£ increases accordingly (see Kallman & Bautista 2001), but the observed equivalent widths of the 
SiXIV line is not much different. This result is explained by the presence of the high-ionized warm absorber at 
log£ > 3, which produces a significant amount of SiXIV but not SiXIII. Consequently, observed variation of the 
Mg and Si absorption line equivalent-widths, is primarily explained by change of the ionization degree of the 
low-ionized warm absorber according to the flux changes. 

We further simultaneously fit the bright/faint spectra in the 1.0-7.5 keV band with the 3-component model 
(Table 5 and 6) by varying normalization of the direct component and ionization degree of the low-ionized warm 
absorber. In stead of allowing the normalization of the absorbed power-law component to be free, we fixed it to 
the value of a constant total normalization minus the direct power-law normalization, in order to confirm the 
anti-correlation between the direct component and absorbed component as observed in Suzaku (Bottom panels 
in Figures 3 and 4). In Table 5, we show results of allowing the ionization degrees of the low-ionized warm 



12 T. Miyakawa et al. 

Table 3. Results of spectral fitting of Chandra/HETG in 1-7.5 keV with the 3-component model for the average spectrum. 



Interstellar absorption 




N H (10 21 cm- 2 ) 


1.5 (fixed) 


W H 




N H (10 23 cm- 2 ) 


2.4(fixcd) 


log £ 


3.64±0.04 


redshift (1(T 3 ) 




W L 




N H (10 21 cm- 2 ) 


3.7(fixed) 


log £ 


2.04±0.03 


redshift (KT 3 ) 


7.4±0.1 


Ni 




10- 2 ph/s/cm 2 at 1 keV 


1.23±0.01 


N 2 




10- 2 ph/s/cm 2 at 1 keV 


<0.2 


W 2 




N H (10 24 cm" 2 ) 


1.6(fixed) 


log £ 


1.57(fixed) 


redshift (1CT 3 ) 


1 (fixed) 


P 




photon index 


1.76±0.01 


E cut (keV) 


160(fixed) 


7V 3 


iVixO.3 






coslncl 


0.866 (fixed) 


iFe 




line E (keV) 


6.34 (fixed) 


sigma (keV) 


0.018 (fixed) 


norm (10 _5 ph/s/cm 2 ) 


1.2±0.4 


EW (eV) 


25±8 


edge E (keV) 


7.11 (fixed) 


MaxTau 


13+ 010 

u - io -0.09 


line E (keV) 


6 98+ 03 


sigma (keV) 


O.Ol(fixed) 


norm (10 -6 ph/s/cm 2 ) 


n +3.9 
°- u -7.6 


reduced chi-square (d.o.f) 


0.92 (5099) 



f "pexrav" model in xspec is used. 
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Table 4. Results of spectral fitting of Chandra/HETGS spectra in 1.15-1.55 kcV and 1.75 - 2.10 keV for the in- 
tensity sliced spectra when only the power-law normalization and the normalization of absorption lines are varied. 

Bright Faint 
Mg XI absorption line 
line E (keV) 1.3424±0.0004 



sigma (10 3 keV) 


1.4±0.5 


norm (10 _5 ph/s/cm 2 ) 


-1.2±0.3 -1.0±0.2 


EW (eV) 


-1.9±0.5 -2.2±0.4 


Mg XII absorption line 


line E (keV) 


1.4621±0.0004 


sigma (10~ 3 keV) 


1.2±0.6 


norm (10 _5 ph/s/cm 2 ) 


-1.0±0.2 -0.4±0.2 


EW (eV) 


-1.7±0.3 -1.0±0.5 


photon index 


0.95±0.07 


norm (10 _3 ph/s/cm 2 at 1 keV ) 


8.6±0.2 6.0±0.2 


reduced chi-squarc 


1.21 (x 2 /d-oi = 1329/1096) 


Bright Faint 



Si XIII absorption line 



line E (keV) 1.8515±0.0005 
sigma ( 10~ 3 kcV) < 1.7 

norm ( 10~ 5 ph/s/cm 2 ) -1.0±0.2 0.4±0.2 

EW (eV) -2.4±0.5 -1.4±0.7 

Si XIV absorption line 
line E (keV) 2.000±0.002 



sigma ( 10 3 kcV) 


c c+1.9 
°-°-2.1 


norm ( 10 -5 ph/s/cm 2 ) 


-1.5±0.4 -1.1±0.4 


EW (eV) 


-3.8±1.0 -3.9±1.4 


photon index 


0.97±0.15 


norm ( 10 -3 ph/s/cm 2 at 1 keV ) 


7.7±0.8 5.7±0.6 


reduced chi-squarc 


1.17 (x 2 /d.oi = 539/459) 
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Table 5. Results of spectral fitting of Chandra/HETG in 1-7.5 keV with the 3-component model 
for the Bright and Faint spectra varying the ionization parameter of the low-ionized warm absorber. 



Bright Faint 


Interstellar absorption 




N H (10 21 cm- 2 ) 


1.5 (fixed) 


W H 




N h (10 23 cm" 2 ) 


2.4(fixed) 


log £ 


3.64±0.05 3.62 ±0.06 


redshift (1CT 3 ) 


1 09+O.24 


W L 




N H (10 21 cm- 2 ) 


3.7(fixed) 


log £ 


1.95±0.04 1.78±0.04 


redshift (1CT 3 ) 








lO -2 ph/s/cm 2 at 1 keV 


1.47±0.01 1.05 ±0.01 


N 2 




1(T 2 ph/s/cm 2 at 1 keV 


1.5 -Nt 


W 2 




N H (10 24 cm- 2 ) 


1.6 (fixed) 


log£ 


1.57(fixed) 


redshift (1CT 3 ) 


1 (fixed) 


P 




photon index 


1.78±0.04 


E cut (keV) 


160(fixed) 


N 3 


iVix0.3 


i?t 




coslncl 


0.866(fixed) 


iFe 




line E (keV) 


6.34 (fixed) 


sigma (keV) 


0.018 (fixed) 


norm (10 _6 ph/s/cm 2 ) 


7.8±4.3 


EW (cV) 


14±8 18 ±10 


edge E (keV) 


7.11 (fixed) 


MaxTau 


o is +0 - 10 


line E (keV) 


6 98+ ' 03 


sigma (keV) 


O.Ol(fixed) 


norm (10~ 6 ph/s/cm 2 ) 


-5 0+ 3 ' 9 


reduced chi-square (d.o.f) 


0.89 (7505) 



f "pexrav" model in xspec is used. 
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Table 6. Results of spectral fitting of Chandra/HETG in 1-7.5 keV with the 3-component model 
for the Bright and Faint spectra varying the column- density of the low-ionized warm absorber. 



Bright Faint 


Interstellar absorption 




N H (10 21 cm- 2 ) 


1.5 (fixed) 


W H 




N h (10 23 cm" 2 ) 


2.4 (fixed) 


log £ 


3.65±0.05 3.62±0.06 


redshift (1CT 3 ) 


i OO+0.23 
L - iiZ -0.31 


W L 




N H (10 21 cm" 2 ) 


3.4±0.4 3.5±0.4 


log £ 


1.54 (fixed) 


redshift (1CT 3 ) 








10- 2 ph/s/cm 2 at 1 keV 


1.57±0.03 1.11 ±0.02 


N 2 




1(T 2 ph/s/cm 2 at 1 keV 


1.6 -iVi 


W 2 




N H (10 24 cm- 2 ) 


1.6 (fixed) 


log£ 


1.57(fixed) 


redshift (1CT 3 ) 


1 (fixed) 


P 




photon index 


i.83±8:8? 


E cut (keV) 


160(fixed) 


N 3 


iVix0.3 


i?t 




coslncl 


0.866(fixed) 


iFe 




line E (keV) 


6.34 (fixed) 


sigma (keV) 


0.018 (fixed) 


norm (10 _6 ph/s/cm 2 ) 


9.1±4.3 


EW (cV) 


17±8 22 ±10 


edge E (keV) 


7.11 (fixed) 


MaxTau 


is+ 011 


line E (keV) 


6.98183$ 


sigma (keV) 


O.Ol(fixed) 


norm (10~ 6 ph/s/cm 2 ) 


-I0.3t 7 7 3 7 


reduced chi-square (d.o.f) 


0.91 (7505) 



f "pexrav" model in xspec is used. 
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data and folded model 




2 5 
Energy (keV) 

Fig. 5. Spectral fit result for the time-averaged HETGS spectra with the 3-componcnt model. 

absorber to be free while the column density is fixed; we see the parameter is significantly variable (log£ = 
1.95±0.04 and 1.78±0.04). Also, ionization degree of the high-ionized warm absorber is made a free parameter, 
but these values are hardly variable (log£ = 3.64±0.05 and 3.62±0.06). 

Because of the presence of constant high-ionized warm absorber, varying column density of the low- ionized 
warm absorber effectively changes the average ionization states. Therefore, in stead of varying ionization degree 
of the low-ionized warm absorber, the bright and faint spectral variation can be equally explained by varying 
its column density. This result is shown in Table 6. 

In summary, Chandra/HETGS data confirm that the spectral variation is described by change of only two 
spectral parameters in the framework of the 3-component model (Equation 2), normalization of the direct 
component (iVi) and either ionization degree or column-density of the low-ionized warm absorber(Wi). 

3.5. Time-variation of the spectral parameters and introduction of the "covering fraction" 

Next, we create a time-sequence of the energy spectra from Suzaku data, and apply the 3-component model. 
We have made energy spectra every 20 ksec, and performed spectral fits to all the sequential 20 ksec spectra 
with the 3 component model, to investigate time variation of the spectral parameters. We chose a 20 ksec 
interval, since variation amplitudes in the 1-10 keV band is maximum at around 20-80 ksec (Inoue et al. 2011). 

In the spectral fits, Ni, N% and Wl are treated as time- variable in Equation (2). We have seen that observed 
variation of Wl is equally represented with change of the column density Nh,l or ionization parameter, £ 
(section 3.3). Here, we vary only Nh.l, while £ is assumed to be invariable throughout the observation. As we 
shall see below, variation of Njj.l is considered more physically reasonable than that of £. The fits are made 
for the three observation sequences separately (obsID:700007010, 700007020 and 700007030). 

Figure 7 shows a correlation between N\ and N2, and an inverse correlation is clearly seen between them, 
as well as in bottoms panels in Figure 3 and 4. We find the slope of the anti-correlation is about unity, which 
suggests that sum of the two parameters is almost constant. Therefore, we introduced the "total normalization" , 
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Mg XI, Mg XII 
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Fig. 6. (Top) Chandra HETGS simultaneous fit in the 1.15—1.55 keV band for the bright spectrum and faint spectrum with 
a power-law plus two negative gaussian model. Mg absorption lines are at 1.34 keV and 1.47 keV. (Bottom) Simultaneous 
fit in the 1.75—2.10 keV band with a power-law plus two negative gaussian model. Si absorption lines are at 1.85 keV and 
2.00 keV. 
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N, as, 
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N = Ni +N 2 . 
Similarly, we introduce the "covering fraction" , a, as 



(3) 
(4) 
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Fig. 7. Correlation between normalizations of the direct component and the absorbed component with the fitting with every 
20 ksec. 



Obviously, 

1 — a = — , 

N ' 



(5) 



thus, instead of (N\,N2,Nh,l), we can equally use (N,a,Nn,L) to describe the observed spectral variation. In 
Figure 8, we show correlations among these three parameters. We can see a good "orthogonality" between the 
two parameters N and a; the total power-law normalization N is not significantly variable, while the covering 
fraction, a, and column-density of the low- ionized absorber, Nh,l, indicate a clear correlation. Namely, spectral 
variation of MCG-6-30-15 at a timescale of 20 ksec is primarily described by a single parameter, either a or 
Nh,l- We remind, however, there is a hint of weak anti-correlation between a and N (top-panel of Figure 8). 
Further sophistication of the model may explain this correlation (section 4.4.4). 

Figure 9 indicates another representation of the variation of spectral parameters. As a function of time, the 
three parameters, N,cx,Nh.l, respectively normalized by their average values, are indicated. We can see N is 
not significantly variable, while a and Nh,l are in sync and much more significantly variable. In the following, 
we will discuss physical meaning of these spectral variations. 
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Fig. 8. Correlation among the covering fraction, a, total normalization, N, and column density of the low-ionized warm 
absorber, Nh,l for the spectral fitting for every 20 kscc. 



4. Discussion 



4-1- Variable partial covering Model 

We have seen that Suzaku and Chandra energy spectra of MCG-6-30-15 are successfully described by the 3- 
component model introduced with Equation (2). The spectral variation indicates that the "total normalization" 
(Equation 3) is not so variable, while the covering fraction (Equation 4) and column-density of the low-ionized 
warm absorber are in sync and more significantly variable (Figure 9). 

How should we interpret this characteristic spectral variation? When we introduced the 3-components in 
Equation (2), we have left the interpretation of the term N2W2 either ionized disk reflection component or 
heavily absorbed component. Meanwhile, we found N2 is clearly anti-correlated with normalization of the 
direct component, N\ (bottom panels of Figures 3 and 4, and Figure 7). On one hand, it is difficult to imagine 
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Fig. 9. Fractional variations of the total normalization, N, column density of the lower-ionized warm absorber, Nh.l, and 
the covering fraction, a, for every 20 ksec. 
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a situation that the disk reflected flux inversely correlates with the directly observed flux. On the other hand, 
if the N2W2 term is due to absorption of the flux emitted from the central X-ray source, the anti-correlation is 
naturally understood if only the partial covering fraction of the central X-ray source, a, is variable, while the 
total flux, N, is constant. The "Difference Variation Function analysis" introduced by Inoue et al. (2011) also 
confirms that the (N, a) set is more orthogonal and fundamental than the (Ni, N2) set. 

Using the total normalization N and the covering fraction a, our 3-component model (2) is rewritten as 
follows: 

F = W H W L (l-a + aW 2 )NP + N 3 RP + I Fe , (6) 

where Wl = cxp — (er(i?, £) Nh.l)- Wc found the column density of the low-ionized absorber Nh,l and a are 
correlated (bottom panel of Figure 8), which is understood as follows: The low- ionized absorber represented by 
Wl is optically thin (Table 1), so that Wl ~ 1 — cr(E,£) Nh.l- On the other hand, the same absorption effect 
can be expressed by a partial covering model with a fixed column-density as, 

W L = 1 - a + a cxp(-a( J B, N^ ed) ) (7) 

«l-aa(£,0<r d) - (8) 
Namely, 

N HiL = aN^ ed \ (9) 

that is seen in the bottom panel of Figure 8. 

4--2. Explanation of the spectral variation 

As we have seen, the observed spectral variation of MCG-6-30-15 is primarily explained by variation of the 
partial covering fraction, a, through equations (6) and (8). Below, we designate our model as the "variable 
partial covering model" . 

In fact, in the framework of the variable partial covering model, the Suzaku intensity-sliced spectra in 1 - 40 
keV is explained mostly by change of the covering fraction with minor change of the total normalisation and 
all the other parameters being invariable. In Figure 10 and Table 7, we show an example of the fit of Suzaku 
intensity-sliced spectra (the same dataset as used in section 3.3). The eight intensity-sliced spectra, besides 
the brightest one, are assumed to have the same total normalization, N. They are fitted simultaneously and 
successfully with the model expressed with Equations (6) and (8) only varying a from 0.63 (dimmest) to null 
(brightest). 

4-3. Explanation of the small variability in the iron energy band 

The well-known small variability in the iron energy band (Section 1) is also nicely explained in the framework 
of the variable partial covering model. In order to see energy dependence of the variation at a specific timcscale, 
we took the following method to calculate the "Difference Variation Function", DVF(E) (Inoue et al. 2011): 
(1) We choose a time-interval of interest, t, to investigate for spectral variations; here, we take 40 ksec. (2) 
Energy spectra are calculated for every t/2, namely, 20 ksec. (3) For every two contiguous energy spectra, we 
recognize the "brighter spectrum" and the "fainter spectrum", B(E) and F(E), respectively. (4) We average all 
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the brighter spectra and all the fainter spectra. (5) Thus, for a given time-period t, we create a single "averaged 
bright spectrum" < B(E) > and a single "averaged faint spectrum" < F(E) >. (6) We calculate difference of 
the averaged bright spectrum and the averaged faint spectrum relative to the average, as 



DVF(E)^ <B ^ > - <F ^ > 



(10) 



<B{E)> + <F{E)>' 

In Figure 11, we show the Difference Variation Function thus calculated at a timescale of 40 ksec. Together, 
we show a model prediction of the Difference Variation Function calculated in the framework of the variable 
partial covering model (Eq. 29 in Inoue et al. 2011), where the covering fraction a is assumed to be the major 
variable parameter, while the total normalization N adds minor fluctuation. Agreement of the observation and 
model calculation is obvious. 

4-. 4- Parameters of the warm absorbers 

Next, we estimate physical parameters of the warm absorbers in Equation (6), Wjj, Wl and W2, which 
are high-ionized warm absorber, thin low-ionized warm absorber and thick partial absorber, respectively. We 
provide suffixes to distinguish column densities and ionization parameters for different warm absorbers. In the 
following, r is distance from the central X-ray source to the absorber, h is a representative thickness of the 
absorbing region along the line of sight, and n is density. 

data and folded model 



0.1 r 
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10- 3 r 




5 - 
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mm ttww 
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Energy (keV) 
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Fig. 10. Suzaku intensity-sliced spectra simultaneously fitted with varying only the covering fraction and all the other 
parameters linked, except the brightest spectrum where the total normalization is also varied. 



First, in general, we note the following relations hold: 
£ = L/nr 2 = ^ ( — ) , where N H = nh, 



NhT \ r 

L fh\ L 

< 



(11) 



(12) 



For the luminosity, we take the representative value 10 43 erg/s, which is valid for order estimation, since the 
flux variation is much less than an order of magnitude. 
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Table 7. Results of spectral fitting in 1—40 keV for the intensity-sliced spectra when 
only the covering fraction is varied (from slice 1 to 8, the dimmest to the brightest). 



state slice 1 slice 2 slice 3 slice 4 slice 5 slice 6 slice 7 slice S 



Interstellar absorption 




Nh (10 21 cm~ 2 ) 


1.9±0.1 






?vr / 1 23 ■ — 2 \ 

Nh (10 cm ) 


2.4 (fixed) 


log i 


3.43±0.03 


Covering Fraction 




a 0.632 


0.356 0.269 0.195 0.181 0.166 0.122 < 0.1 


W L 




ivr( fixed) /t^2^ — 2\ 

N H L (10 cm ■ a ) 




Nh 




log 5 


1.88±0.03 






10 pn/s/cm at 1 kev 


1.80 2.71 


Ni 


(1 — a) N 


N 2 


a N 


W 2 




Nh (10 cm ) 


1.6 (fixed) 


log £ 


1.63±0.05 


P 




photon index 


1.88±0.01 


E cui (kcV) 


160 (fixed) 






pexrav K (10 — 3 ) 


0.3 x TVi (for slice5) 


h" 




coslncl 


0.866 (fixed) 






line E (keV) 


6.35 (fixed) 


sigma (keV) 


0.01 (fixed) 


/in — 5 1 / / 2n 

norm (10 pn/s/cm J 


1.4±0.2 


edge E (keV) 


7.11 (fixed) 


MaxTau 


0.08±0.01 


line E (keV) 


7.0 (fixed) 


sigma (keV) 


0.01 (fixed) 


norm (10~ 6 ph/s/cm 2 ) 


-4.8±1.7 


line E (keV) 


2.35±0.01 


sigma (keV) 


0.01 (fixed) 


norm (10 _5 ph/s/cm 2 ) 


-1.9±0.3 


Reduced chisq (d.o.f) 


0.49 (1149) 



f "pexrav" model in xspec is used. 
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Fig. 11. Energy dependence of the spectral variation ("Difference Variation Function"; blue) and model simulation (red) in 
the framework of the variable partial covering model. This figure is the same one as Figure 4 in Inoue et al. (2011). 

.4--1- High-ionized warm absorber, Wh 

From model fitting (Table 7), £ H ^ 10 3 ' 4 erg cm/s, N H . H ~ 10 234 cnr 2 . Hence, 



L 



H 



< 



L 



10 



16.2 



cm. (13) 

Nh.h^h \th J Nh,h^h 
Our analysis has indicated that the high- ionized warm absorber is constant, while partial covering fraction is 
variable. This suggests the variable partial absorbers be embedded in the largely extended, static high-ionized 
warm absorber. Hence, we assume hn ~ rjj, and take ~ 10 16 cm and nn ~ 10 7 cm -3 for typical parameters 
of the high-ionized warm absorbers. 
4-4-%- Thick partial absorber, W 2 

From model fitting (Table 7), £> ^ 10 16 erg cm/s, N H ,2 ^ 10 24 2 cm- 2 . Hence, 



L 



r-2 



h 2 



< 



L 



10 



17.2 



(14) 



Here, we may further constrain the parameters of the partial absorber from variation timescale of the partially 
absorbed component. In our variable partial covering model, typical variation timescale (~ 10 5 sec; Figure 9) 
corresponds to the crossing time of the X-ray absorbing clouds in front of the X-ray source. We will call these 
clouds as the "low-ionized clouds" , hereafter. Assuming that size of the X-ray source and that of a low-ionized 
cloud have similar dimensions, the crossing timescale is expressed as ~ h-2.IV , where V is the velocity of a 
low-ionized cloud. Hence, 

V 



h 2 ~ V x 10° ~ 10 1 



10 9 cm/s 

where V is normalized to a typical velocity of the broad line region (BLR) clouds. Combining h 2 and Nh,2 
10 24 cm -2 , we estimate parameters of the partial absorbers as follows: 



(15) 



L 



l'2 



1/2 



10 157 cm, 



(16) 



25 



n 2 ~ 10 1U cm 



(17) 



Note that a low-ionized cloud may have internal ionization structures, which we will sec below. 
4-4-3- Thin low-ionized warm absorber, Wl 

In our variable partial covering model, the low- ionized warm absorber (Wl) and the partial absorber (W2) are 
closely associated. Presumably, they are different parts of the low-ionized clouds. We point out that presence 
of such internal strcutures in the BLR cluds is also suggested in Mrk 766 (Risaliti et al. 2011). 

In this model, the low-ionized warm absorber has a fixed column-density jvj^" ed ' ) with a variable partial 
covering factor (Equation 9); we simply write this column density as Nh,l below. 

From model fitting (Table 7), N h ,l ~ 10 22 1 cm" 2 and £l ~ 10 19 erg cm/s. 

If we take the distance of the low- ionized clouds r 2 , 



Note that values of ul and £l are between those of the high-ionized absorber (n# ~ 10 7 cm -3 and ~ 10 3 4 ) 
and the partial absorber (n 2 ~ 10 10 cm -2 and £2 ~ 10 ), respectively The thickness Hl is much smaller 
than that of the partial absorber h 2 ~ 10 14 cm. Therefore, the low-ionized absorber is considered to be in the 
boundary layer between the low-ionized cloud and the high-ionized warm absorber. We may call this putative 
layer as the "cloud envelope" . 

Meanwhile, it is known that there takes place thermal instability in the range of £ = 10 2 ~ 10 3 erg cm/s 
(Reynolds & Fabian 1995). Therefore, it is reasonable that the cloud envelope continuously connects the partial 
absorber in the low-ionized cloud and the surrounding constant high-ionized warm absorber. 
4-4-4- Putative X-ray blocker 

We notice the Thomson optical depth (~ Nh /10 24 2 ) of the partial absorber is about unity. Presumably, 
this is not just a coincidence, but suggests that the low-ionized cloud has internal density gradients extending 
up to more than unit Thomson optical depth, so that attenuation with continuously different optical depths is 
approximated by unit Thomson optical depth. Thus, the low-ionized cloud may well have central, cold Thomson 
thick cores (Nh > 10 24 2 cm -2 ), which are completely opaque to incoming X-rays. We may not see these "X-ray 
blockers" directly, but if we take into account the effect of the X-ray blockers, the weak anti-correlation between 
the total normalization N and the covering fraction a (top-panel of Figure 8) might be explained. 

4-5. Origin of the Variable Partial Covering 

Finally, we present a physical picture surrounding the X-ray source to explain the observations. In Table 8, 
we summarize parameters of the absorbers in our model. Basic points of our variable partial covering model 
are the following: 

1. Observed X-ray flux and spectral variability is primarily caused by random passages of the low- ionized 
clouds in front of the central X-ray source. A typical size of each low-ionized cloud is h 2 ~ 10 14 cm and a 
typical velocity is V ~ 10 9 cm/s (Equation 15). 



~ 10 97 cm- 3 , 



(18) 



h L = N H , L /n L ~ 10 12 - 4 cm. 



(19) 
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Table 8. Summary of Parameters of the Absorbers 



Absorbers r [cm] log Nu log£ h [cm] 



High-ionized warm absorber, Wh 

Envelope Thin low-ionized warm absorber, Wl 
Low-ionized cloud ^ Main body Thick partial absorber, W 2 

Core X-ray Blocker 



10 16 


23.4 


3.4 


10 16 


10 7 




f 22.1 


1.9 


lfJ 12.4 


1Q 9.7 


10 15.7 


24.2 


1.6 


10 14 


10 10 




>24.2 


< 1.6 


«10 14 


>10 10 



2. Central X-ray source is not always fully blocked, nor fully exposed. The covering fraction is variable 
between null to ~ 0.6 (Table 7), and it fluctuates by an amplitude of ~ 50 % (Figure 9). 

We may estimate the average covering fraction C ave from intensity-sliced spectra as C ave = > Cj x U/T, 

— J i— 1 

where ti and T are exposure time for each spectrum and the total observation, respectively. We thereby obtain 
Cave ~ 0.3. According to Blandford et al. (1991), the integrated covering fraction is about 0.1 for the broad 
line region (BLR). Our results suggest that the partial absorption takes place in the BLR, and the low- ionized 
clouds corresponds to the BLR blobs. 

We may estimate number of the absorbing clouds in the field of view, N c i oua -, and size of the central X-ray 
source, r x , as follows. Assuming that fluctuation of the covering fraction, w 0.5, is determined by statistical 
fluctuation of the number of clouds, AN c i ou d/N c i oua ; N c i ou d is estimated as w 4. The average covering fraction, 
rj 0.3, may be written as N c i ou dir(h2/2) 2 /Trr x 2 , thus, r x « 2/i2 w2x 10 14 cm. 

Typically, BLR blobs have a velocity of V ~ 10 9 cm/sec. Assuming the Kepler velocity has a similar value, 
r = (c? /2V 2 )r s ~ 500 r s , where r s is the Schwarzschild radius. Identifying this with r2 ~ 10 15 7 cm (Equation 16), 
the central black hole mass of MCG-6-30-15 is ~ 3 x 10 7 M Q . Thus, the central X-ray source (r x ~ 2 x 10 14 cm) 
is extended to ~ 20 r s around the black hole. Besides variation of the BLR coluds, the X-ray source is variable 
on a timescale down to ~ 1 ksec (Inoue et al. 2011), which may be compared with the free fall time, r/yj ~@M.. 
Then, the intrinsice source variation arises at a few times the Schwarzshild radius, around the innermost stable 
circular orbit. 

The low-ionized clouds are presumably located along the equatorial plane of the accretion disk, assuming 
they are related to the accretion flow. Thus we are supposedly observing MCG-6-30-15 close to cdgc-on.A 
similar partial covering model with two layers of absorbers is proposed for the Galactic black hole candidate 
GRO J1655-40 to explain spectral varaitions during the dips (Tanaka, Ueda and Boiler 2003). Solid angle of 
the low-ionized clouds seen from the X-ray source should be much smaller than 47r, otherwise we will have to 
see strong fluorescent lines from the clouds. If a thick absorber with Nh ^ 10 24 cm" 2 covers the power-law 
component completely (51 ~ 47r), equivalent width of the iron line to the reflection component is ~ 1 kcV. 
Considering C ave ~ 0.3 and the ratio of the reflection component to the direct component, 10 to 1, the expected 
iron line equivalent width is ~ 30 eV. As mentioned in subsection 5.5.2, we found the fit improves if a mildly 
broad disk-line is added, where the best-fit inner radius and the equivalent width are > 200 r g and 40 ± 20 eV, 
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respectively (Table 2). These mildly broadened, weak iron emission lines may well be expected from fluorescence 
in the BLR clouds. 

Figure 12 gives a schematic view of internal structure of the low-ionized cloud and the variable partial covering 
model surrounding the black hole. The low-ionized cloud (BLR cloud) has internal ionization structure. A 
putative blocker in the core would block the incoming X-rays completely. The main body of the low-ionized 
cloud is the thick partial absorber surrounding the core. Still the outer cloud envelope corresponds to the 
thin low-ionized warm absorber. Central X-ray emission region around the black hole is mildly extended, and 
partially covered by these low-ionized clouds. Variation of the covering fraction explains most of the observed 
flux and spectral variations. 




Fig. 12. Schematic picture of the variable partial covering model for MCG-6-30-15 and structure of the absorbing cloud. 

5. Conclusion 

We have analyzed Suzaku and Chandra archival data of MCG-6-30-15 and constructed a new model to explain 
the observed spectral variation. Our main conclusions are summarized below: 

1. We have shown that the observed energy spectra and spectral variation can be explained by the "3- 
component model" (Equation 2) that includes (1) a direct power-law component, (2) a heavily absorbed 
power-law component by thich photoionized material, and (3) a cold disk reflection component far from 
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the black hole. The first two components are affected by two warm abosrbcrs having different ionization 
states. 

2. The ionized iron K-edge of the heavily absorbed component and spectral curvature due to the warm ab- 
sorbers explain most of the seemingly broad "disk line" spectral feature. Consequently, general relativistic 
interpretation of the broad iron emission line feature in MCG-6-30-15 is not confirmed. The claim that 
MCG-6-30-15 is a Kerr black hole with nearly extreme rotation assuming the relativistic line-broadening 
(e.g., Miniutti et al. 2007) is questionable. 

3. We propose the "variable partial covering model", in which a central, moderately extended X-ray source 
is partially covered by variable low-ionized clouds in the line of sight. In this model, observed X-ray 
spectral variation is primarily caused by change of the partial covering fraction of the central source. 
These absorbing clouds presumably correspond to the fast-moving broad line region (BLR) clouds. 
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